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Reversing the magnetization of a ferromagnet by spin transfer from a current, rather than by 
applying a magnetic field, is the central idea of an extensive current research. After a review 
of our experiments of current-induced magnetization reversal in Co/Cu/Co trilayered pillars, we 
present the model we have worked out for the calculation of the current-induced torque and the 
interpretation of the experiments. 

PACS numbers: 



The concept of magnetization reversal by spin trans- 
fer from a spin-polarized current was introduced in 
1996 by Slonczewski [1]. Similar ideas of spin trans- 
fer had also appeared in the earlier work of Berger [2] 
on current-induced domain wall motion. Convincing ex- 
periments of magnetization reversal by spin transfer on 
pillar-shaped multilayers [3-6] , nanowires [7] or nanocon- 
tacts [8] have been recently performed and several the- 
oretical approaches, extending the initial theory, have 
also been developed [9-19]. From the application point 
of view, magnetization reversal by spin transfer can be 
of great interest to switch spintronic devices (MRAM 
for example), especially if the required current density 
- presently around 10'' A/cm^ - can be reduced by ap- 
proximately an order of magnitude. 

We present a summary of our experiments on 
Co/Cu/Co pillars, describe a calculation model for the 
critical currents as a function of - mainly - CPP-GMR 
data and we discuss its application to experiments. 

I. Experiments 

We present experiments on pillar-shaped Col (2.5 
nm)/Cu(10 nm)/Co2(15 nm) trilayers. The submicronic 
(200 X 600 nm^) pillars are fabricated by e-beam lithog- 
raphy [5]. The CCP-GMR of the trilayer is used to detect 
the changes of the magnetic configuration (the difference 
between the resistances of the P and AP configurations 
is about 1 mfi). For all the experiments we describe, 
the initial magnetic configuration is a parallel (P) one, 
with the magnetic moments of the Co layers along the 
positive direction of an axis parallel to the long side of 
the rectangular pillar. A field ffappi is applied along the 
positive direction of this axis (thus stabilizing this ini- 
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FIG. 1: Resistance vs. dc current: (a) sample 1 for Happi = 
(black) and Happi = 125 Oe (grey); (b) sample 2 for ffappi = 
(black), ffappi = + 500 Oe (grey) and Happi = + 5000 Oe 
(dotted line). 



tial P magnetic configuration). We record the variation 
of the resistance [R) as the current / is increased or de- 
creased (positive / means electrons going from the thick 
Co layer to the thin one). The results we report here 
are obtained at 30 K (the critical currents are smaller at 
room temperature). 

In Fig. 1(a), we present a typical variation of the 
resistance i? as a function of the current, for iJappi = 
and -f 125 Oe. Starting from a P configuration at 
/ = and increasing the current to positive values, we 
observe only a small progressive and reversible increase 
of the resistance, which can be ascribed to Joule heating 
(this has also been seen in all other experiments on pil- 
lars [3-6] when the current density reaches the range of 
10^ A/cm^). In contrast, when the current is negative 
and at a critical value Iq^^^, there is an irreversible 
jump of the resistance (Ai? « 1 mil), which corresponds 
to a transition from the P to the AP configuration (re- 
versal of the magnetic moment of the thin Co layer). 
The trilayer then remains in this high resistance state 
(the Rap {I) curve) until the current is reversed and in- 
creased to the critical value 1^^^^ , where the resistance 
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FIG. 2: Instability lines of the P and AP configurations (sam- 
ple 1). The P configuration is stable above line 1 and unsta- 
ble below. The AP one is stable below line 2 and unstable 
above. At low field [regime A), the stability zones of P (blue) 
and AP (yellow) overlap between Imes 1 and 2 (stripes). At 
high field [regime B), there is a zone (green) between lines 1 
and 2 where none of the P and AP configurations is stable. 
Equations of lines 1 and 2 are derived from a LLG equation 
for uniaxial anisotropy If an [18]. The magnetic field includes 
f^appi. and, possibly, interlayer coupling fields, lines 1 and 2 
cross at about i/an. 



drops back to the Rp{I) curve. This type of hysteretic 
R{I) cycle is characteristic of the magnetization reversal 
by spin injection in regime A. 

For iJappi — 0, I^^^^ = — 15 mA (current density 
jP^AP _ X 10^ A/cm2) and 1^^^^ ^ -h 14 mA 
i.Jc^'^^ - + 1-17 X 10'' A/cm^). A positive field, which 
stabilizes the P configuration, shifts slightly the critical 
currents; \Ic^^^\ increases and I^^^^ decreases (note 
that the relatively larger shift of I^^^^ at 125 Oe in 
Fig. 1(a) is specific to the approach to the crossover to 
regime B at about 150 Oe). 

The R{I) curve for i/appi = + 500 Oe, shown in Fig. 
1(b), illustrates the different behavior when the applied 
field is higher [regime B). Starting from / = in a P 
configuration (on the Rp{I) curve), a large enough neg- 
ative current still induces a transition from P to AP, but 
now this transition is progressive and reversible. The 
R[I) curve departs from the Rp{I) curve at ifuvi-t^ — 

— 25 mA [jftZf^ = - 2.08 x 10^ A/cm^) and catches up 
the Rap{I) curve only at 4^7^^ ^ - 45 mA (j^^^^ = 

— 3.75 x 10^ A/cm^). On the way back, reversibly, R[I) 
departs from Rap{1) at ifta^^ = ler^^^ ^ - 45 mA 
and reaches finally Rp[I) at ltn<r = ^ ^ft^rt^ = 

— 25 mA. At higher field, the transition is similarly pro- 
gressive and reversible, but occurs in a higher negative 
current range. Finally, for very large applied field (-ffappi 
= 5000 Oe), the transition is out of our experimental 
current range, and the recorded curve is simply Rp[I). 

The experimental results presented above can be sum- 
marized by the diagram of Fig. 2. This type of diagram is 
obtained [18] by introducing the current-induced torque 
into a Landau-Lifshitz-Gilbert (LLG) motion equation to 
study the stability /instability of the moment of the mag- 
netic thin layer (the moment of the thick layer supposed 



being pinned). The P configuration is stable above line 
1 and unstable below. The AP configuration is stable 
below line 2 and unstable above. 

Regime A corresponds to Hg^ppi smaller than the field 
at which line 2 crosses line 1. In this regime, there is an 
overlap between the stability regions of P and AP. Start- 
ing from a P configuration at zero current and moving 
downward on a vertical line, the P configuration becomes 
unstable at the negative current I^^^^ corresponding 
to the crossing point with line 1. As this point in the 
stability region of the AP configuration, the unstable P 
configuration can switch directly to the stable AP config- 
uration. On the way back, the AP configuration remains 
stable until the crossing point with line 2 at I^^^^ (pos- 
itive), where it can switch directly to a stable P con- 
figuration. This accounts for the direct transitions and 
hysteretic behavior of regime A in Fig. 1(a). 

In regime B, for i?appi above the crossing point of lines 
1 and 2, none of the P and AP configurations is stable 
in the region between lines 1 and 2. Going down along a 
vertical line, the P configuration becomes unstable at the 
crossing point with line 1 [iftart'^) and the system de- 
parts from this configuration. But the AP configuration 
is still unstable at this current and can be reached only 
at the crossing point with line 2 {ler^^^)- the way 
back, reversibly, the AP configuration becomes unstable 
at the crossing point with line 2 [I^art^ = lend^^)^ but 
a stable P configuration is reached only at the crossing 
point with line 1 {ifj^^^ — I^tarf^)- This accounts for 
the behavior of Fig. 1(b). The state of the system dur- 
ing the progressive transition between P and AP can be 
described as a state of maintained precession. 

The critical lines of the diagram of Fig. 2 can also be 
derived from the variation of R along a horizontal line, 
for example from the i?(i?appi) curves of Fig. 3 for sample 
2. The i?(i?appi) curve for / ~ +50 mA is flat, i.e. there 
is no GMR. This is because, along an horizontal line in 
the upper part of the diagram of Fig. 2, the P configura- 
tion is always stable. For negative current, on the other 
hand, the i?(-ffappi) curves mimic the GMR curves of an 
antiferromagnetically coupled trilayer, in which the anti- 
ferromagnetic coupling would increase when the current 
becomes more negative. This can be expected from the 
diagram of Fig. 2. For example, starting from high field 
at / = — 50 niA, the upturn from the baseline at about 
Happi — + 5600 Oe indicates the beginning of the pro- 
gressive transition from P to AP at the crossing point 
with line 1. As -ffappi is decreased further, the progres- 
sive (and reversible) increase of R reficcts the progressive 
crossover from P to AP on a horizontal line between line 1 
at 5600 Oe and line 2 at a. field in the range 100—200 Oe. 
When the moment of the thick Co layer is reversed in a 
small negative field, the P configuration being unstable 
and the AP stable in this region of the diagram, the mo- 
ment of the thin layer is also reversed to restore the AP 
configuration, so that R is practically not affected by the 
coupled reversal of both layers. 

We conclude that the main features of the experimental 
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FIG. 3: Resistance vs. applied magnetic field in sample 2 
for / = — 50 mA, — 40 mA, and — 30 mA. For clarity, the 
curves have been shifted vertically to have the same high field 
baseline, inset ■.R{H) for I = + 50 mA. 

results fit into tlie frame of tlie diagram of Fig. 2. In Sec- 
tion IV, we discus more quantitatively the influence of pa- 
rameters such as layer thicknesses, spin diffusion length, 
etc. The final remark of this Section is that the phase 
diagram of Fig. 2 comes from an oversimplified model 
assuming that the only current-induced excitations are 
precessions of a global magnetization vector due to trans- 
verse spin transfer. Several types of additional effects can 
be expected from non-uniform precessions, or, more gen- 
erally, from other modes of current-induced excitations. 
For example, excitation of magnons is probably a signif- 
icant dissipation mechanism in the stage of maintained 
precession and also a dissipation channel of the longitu- 
dinal spin accumulation at high current density. These 
additional excitation modes should also be reflected in 
the resistivity and probably account for some features of 
the experimental results [3-8] which arc not described by 
the scheme of Fig. 2 for pure rotations. Others effects [6] 
are also expected from exchange or dipolar interlayer cou- 
plings which can play the same role as the applied field 
in Fig.2. 

II. Theoretical Model 

The magnetization of a magnetic layer can be reversed 
by spin transfer if the spin polarization of the injected 
current and the magnetization of the layer are non- 
colinear. In a multilayered structure this requires a non- 
colinear configuration of the magnetizations of the differ- 
ent layers. The transfer from an obliquely polarized spin 
current running into a magnetic layer is associated with 
the alignment of the polarization of the current inside 
the layer along the magnetization axis. If the current- 
layer interaction is spin conserving (exchange-like), this 
implies that the transverse component of the spin cur- 
rent is absorbed and transferred to the layer. This is the 
spin transfer concept introduced by Slonczewski [1] . The 
contribution of this transfer to the motion equation of 
the total spin S of the layer is written as: 

{dS/dt)j = absorbed transv. spin current (1) 



or, in other words, a torque equal to the absorbed spin 
current multiplied by h is acting on the magnetic moment 
of the layer. 

Several mechanisms contribute to the transfer of the 
transverse component of a spin current running into a 
magnetic layer [12]. First, due to the spin dependence of 
the reflection/transmission process at the interface with 
a ferromagnet, the transverse component is reduced and 
rotated in the transmitted spin current. What remains 
of transverse component then disappears (is transferred) 
by incoherent precession of the electron spins in the ex- 
change field of the ferromagnet. Ab-initio calculations 
[12] show that, for a metal like Co, the transverse spin 
current is almost completely absorbed at a distance of the 
order of 1 nm from the interface. In these conditions, the 
spin transfer is a quasi-interfacial effect and, in our calcu- 
lation, is expressed by interface boundary conditions (in 
the same way as interface resistances are introduced in 
boundary conditions for the theory of CPP-GMR [20]). 
On the other hand, the longitudinal component of the 
spin current in the magnetic layers and all its compo- 
nents in the nonmagnetic layers vary at the much longer 
scale of the spin diffusion length l^i (60 nm in Co, about 
1 /im in Cu). They can be calculated by solving diffu- 
sive transport equations for the entire structure, as in the 
theory of the CPP-GMR. An essential point is that, for 
a non-colinear configuration with different orientations 
of the longitudinal axes in different layers, the longitu- 
dinal and transverse components of the spin current are 
inter-twined from one layer to the next one, so that a 
global solution for both the longitudinal and transverse 
component and for the entire structure is required. 

The calculation of our model can be summarized as fol- 
lows. We consider a Ni/Fi/N/F2/N7j structure, where 
Fi (thin) and F2 (thick) are ferromagnetic layers sepa- 
rated by a thick nonmagnetic layer N. and are 
two semi-infinite nonmagnetic layers (leads). For sim- 
plicity we assume that Fi and F2 (N, and N/j) are 
made of the same ferromagnetic (nonmagnetic) material. 
The current is along the x axis perpendicular to the lay- 
ers. m(x) and j{x) are the 2x2 matrices representing 
respectively the spin accumulation and the current den- 
sity: 

= '^T+ j,n,x{x)(Tx + jrn,yix)dy + jm,z{x)az 

m{x) = mx{x)dx + my{x)dy + mz{x)dz (2) 

where CTx, and dz are the three Pauli matrices and 
I is the unitary matrix. Spin accumulation and current 
are defined as in Ref [13]. If we call Zi the local spin po- 
larization axis {zi = zi in Fi, Zi = Z2 in F2), m^. (jm.zj 
is the longitudinal component of the spin accumulation 
vector m (spin current vector jm), nix^ and my^ {jm.xi 
and jm,vi) are the transverse components of m {jm)- 

To derive the critical currents for the instability of the 
P and AP configurations, we need only to calculate the 
current-induced torque in the simple limit where the an- 
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gle between the magnetizations of the magnetic layers is 
small or close to tt (0 or tt — 9, with 9 small). The first 
step, before introducing the small angle 9, is the calcula- 
tion of the longitudinal spin current jmz and spin accu- 
mulation niz in a colinear configuration [9 — Q). This is 
done by using the standard diffusive transport equations 
of the theory of the CPP-GMR with parameters (spin de- 
pendent interface resistances, interface spin memory loss 
coefficient, spin diffusion lengths, etc) derived from CPP- 
GMR experiments [21,22]. An example of the result for 
the P configuration of a Co/Cu/Co trilayer is shown at 
the top left of Fig. 4. In the bottom part of Fig. 4, we 
represent the situation when a small deviation 9 from the 
parallel colinear configuration above is introduced. The 
spin accumulation in the Cu spacer layer is a constant 
vector m (as, generally, ^cu ^ ^S")- With respect to the 
colinear configuration, the amplitude of m has changed 
by a quantity of the first order in 9 (we omit this part 
of the calculation). However, to calculate the torque at 
first order in 9, we can neglect this change and assume 
m| — rriQ^, where is the spin accumulation iriz in 
Cu for the P colinear configuration. On the other hand, 
m cannot be parallel to both zi and Z2, and its orien- 
tation in the frame of the thin layer is characterized by 
the unknown angles 9ra (of the order of 9) and x- These 
angles will be determined later by self-consistency condi- 
tions for the whole structure. The key point, explaining 
the injection of a large transverse spin current into the 
thin magnetic layer, is the discontinuity of transverse spin 
accumulation between the two sides of the interface be- 
tween Cu and Col, |m| = 9m.'m^^ in Cu and |m| = in 
Col. This is equivalent to a huge gradient of spin accu- 
mulation and generates a large transverse spin diffusion 
current running into the interface where it is absorbed 
or reflected. A straightforward angular integration, illus- 
trated at the top right of Fig. 4, gives for the incoming 
transverse spin current: 



1 



(3) 



where j™*! = + ^3™% and Wi? is the Fermi velocity. 
Eq.(|3Jl holds for a spacer thickness of the order of the 
mean free path or larger. A part of this incoming trans- 
verse spin current is reflected into Cu at the Cu/Col 
interface. The remaining part absorbed in the interfacial 
precession zone can be written as = te^^f^^ , where 
the coefficient t and the rotation angle e have been cal- 
culated [12] for a large number of interfaces. This leads 
to: 



■abs. 



(4) 



For thinner spacer layers, a contribution to the diffusion 
current comes also from the thick magnetic layer and 
includes an additional term in m^^ [16]. The scale 
of the transverse spin current of Eq.Q is the product 
m^^vp (or m^^vp around the AP state), where rn^^ is 
controlled by the spin relaxation in the system, rn^^vp is 



of the order of (je/e)(/sf/A), where (laf/X) is a mean value 
of the ratio of the spin diffusion length (SDL) to the mean 
free path (MFP) in the structure (including the leads), 
and can be definitely larger than the charge current je/e- 
In most cases, the transverse spin current of Eq.Q) will 

be larger than the current ^mij^'^M ^ corresponding to the 
projection of the longitudinal spin current in the colinear 
configuration (the diffusion spin current coming from the 
gradient of spin accumulation). 

The unknown angles 9m and x are calculated [16] by 
imposing a global cancellation of the transverse spin cur- 
rents outgoing from or reflected into the spacer layer. In 
the case of a small deviation 9 from the P configuration, 
for example, this leads to 9m = 9/2 and x — ■""Z^, and 
finally, from Eq.l^, to the following general expression 
of the torque T^: 
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xMi A (Ml A Ma) 



(5) 



with a similar expression for (Mi and M2 are 

unit vectors along the magnetizations, ni^^ and j^^ 
are the spin accumulation and current at the Cu/Co2 
interface in the colinear configuration). As ab-initio 
calculations have shown that, for most interfaces be- 
tween classical magnetic and nonmagnetic metals [12], 
t is always close to 1 and e very small (t = 0.92 
and e smaller than 3 x 10^^ for Cu(lll)/Co, for ex- 
ample), we have supposed t = I, e = and kept 
only the term Mi A (Mi A M2) in an expression of 
the form [cos(e)Mi A (Mi A M2) + sin(e)Mi A M2] (as- 
suming e = is equivalent to neglecting the small imag- 
inary parts of the mixing conductances in circuit the- 
ory [15]). In. Eq.lO we have also neglected the interfa- 
cial memory loss of transverse spin by spin-orbit effects 
(the longitudinal spin memory loss at the interfaces [21] 
is already taken into account in the calculation of m and 
im in the colinear configuration). 

The important feature in Eq.(|3Jl is the relation of the 
torque at small angle to the spin accumulation m and 
spin current jm calculated for the P and AP colinear 
configurations. We emphasize that, due to the relevant 
length scale of this calculation, the result for F involves 
the entire structure (including a length of the order of 

P(AP) 

the SDL in the leads). The spin currents jm cu ^^"^ only 
a fraction of the charge current je/e. In contrast the 
terms vpm^^^^\ reflecting the diffusion currents gener- 
ated by the transverse spin accumulation discontinuities 
in a non-colinear system, are of the order of {je/e){lsf/X) 
and can be larger than je/e (a special case, however, 
is that of a P configuration of a symmetric structure. 



for which 



0). Enhancing the spin accumulation 



and increasing its ratio to the current je is certainly the 
most promising way to reduce the critical current, for 
example with materials in which a higher spin accumu- 
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FIG. 4; Top left: Profile of the spin current jm.z 
and spin accumulation niz calculated from diffusive CPP- 
transport equations and CPP-GMR data for a (Cu/Col 2.5 
nm/Cu 10 nm/Co2 oo) structure in a parallel colinear config- 
uration with an electron current (je) going to the left. Bot- 
tom: For a small angle 9 between the polarization axes z\ 
(vertical) and Z2 of the same structure, 3D sketch represent- 
ing the spin accumulation m in the Cu layer (|m| = mcu oi 
the colinear configuration), its transverse component mx and 
the transverse component of the induced spin currents diffus- 
ing to, refiected from and absorbed by the Col layer. The 
angles 9m and Xm characterize the orientation of the vector 
m in the frame of Col. Top right: Schematic illustrating the 
calculation of the transverse spin diffusion current generated 
by the transverse spin accumulation on the Cu side of the 
Col/Cu interface. 



lation splitting can be expected (magnetic semiconduc- 
tors ?). This dependence on SDL and "amplification" is 
also taken into account in the model of Stiles and Zang- 
will [11,12] or Kovalev et al. [15], and in recent calcula- 
tions of Slonczewski [10]. This "amplification" also turns 
out in the model of Shpiro et al. [14] for the opposite 
limit of non-interfacial transfer. The main difference be- 
tween the two limits is the equal importance of the terms 
Ml A (Ml A M2) and Mi A M2 in the torque of Shpiro 
et al. [14]. We will see below that the experimental criti- 
cal line diagram of Fig. 2 indicates a largely predominant 
Ml A (Ml A M2) torque term. 

III. Discussion and Conclusion 

Our expression of the torque, Eq.Q, can be applied 
to the interpretation of the experimental results. 

(a) If the torque of Eq.® is written as G^^'^^'je x 
Ml A (Ml A M2) and, when the excitation can only be 
an uniform precession, the critical currents at zero field 
are expressed as [3,17,18]: 



■P^AP _ 

Jc 

■AP~^P 

Jc 



^70 



{Ha 



2ttM) 
2ttM) 



where a is the Gilbert coefRcient, iJan is the anisotropy 
field and M the magnetization. By using experimental 
data (interface resistances, interface spin memory loss co- 
efRcient, SDL, etc) from CPP-GMR experiments [21,22] 
to calculate the spin accumulation in the Co/Cu/Co 
trilayer and then r^^"*^) and G^^'^'^' from Eq.©, we 
obtain a reasonable agreement with our experiments: 
2.8 X 10^ A/cm^ (exp.: - 1.25 x lO^A/cm^) 



■P^AP _ 

Jc ~ 
and jn^^^ 



1.05 X 10^ A/cm2 (exp.: + 1.17 x lO'^ 



(6) 



What can be also predicted for the critical currents 
is : i) their proportionality to the thickness of the thin 
magnetic layer (this follows from the assumption of in- 
terfacial transfer and has been already observed [3]); ii) 
their decrease as the thickness of the thick magnetic layer 
increases, with saturation at a minimum level when the 
thickness exceeds the SDL (60 nm in Co at low tempera- 
tures, for example); iii) their increase (at the scale of the 
mean free path in the spacer) when the spacer thickness 
increases; iv) their definite dependence on the SDL in the 
layers and leads. 



(b) In finite applied field, a diagram of the type 
of Fig. 2, with a crossover between the two regimes 
around H — H^n, is expected for a torque of the form 
Ml A (Ml A M2). The equations of the critical lines 
and a fit with experimental data has been presented else- 
where [18]. The diagram expected for a torque Mi A M2 
does not include a zone where both the P and AP con- 
figurations are unstable {regime B with progressive and 
reversible transition) and cannot be fitted with the ex- 
periments on Co/Cu/Co trilayers. 



Although the spin transfer effect begins to be better 
understood, the possibility of reducing sufficiently the 
critical currents for practical applications is still a pend- 
ing question. For conventional ferromagnetic metals (Co, 
etc) and from numerical applications of the model of this 
paper [16], some reduction seems possible but probably 
by less than an order of magnitude. As we have pointed 
out, a stronger reduction might be obtained with other 
types of magnetic materials permitting higher spin ac- 
cumulations. On the other hand, another type of spin 
transfer effect is the current-induced domain wall mo- 
tion [2] . According to recent experimental results of do- 
main wall motion with relatively small current densi- 
ties [24] , this should be also a promising way for current- 
induced switching. 
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